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T-H Abstract 

o 

C\| We present a calculation of the leading order hadronic contribution to the 

anomalous magnetic moment of the muon for a dynamical simulation of 2+1 flavour 
QCD using domain wall fermions. The electromagnetic 2-point function is evaluated 
on the RBC-UKQCD lattice gauge configurations and this is fitted to a continuous 
form motivated by models of vector dominance. We determine a robust and reliable 
technique for performing this fit, allowing us to extract the most accurate results 
possible from our ensembles. This combined with data at very light quark masses 

p produces the result 

(D af^""^ = 641(33) (32) x 10"^° 

at the physical point, where the first uncertainty is statistical, and the second is an 
1—1 estimate of systematics, which is in agreement with previous results. We outline 

^ various methods by which this calculation can and will be improved in order to 

compete with the accuracy of alternative techniques of deducing this quantity from 



■^ experimental scattering data. 



o 



1 Introduction 

The anomalous magnetic moment a of a lepton, is half the discrepancy from 2 (a = ^^) 
of g, the gyromagnetic ratio or Lande ^^-factor, which relates the spin S of the lepton to 
its magnetic moment /I as 

^i = 9,^s. (1.1) 

It is given the name "anomalous" because it is a purely quantum effect and so is zero in 
a classical theory. 

The one-loop computation of the electron anomalous magnetic moment a^. by 
Schwinger [T] was one of the first such calculations, and provided strong evidence in 
support of the young theory of quantum electrodynamics (QED) by explaining observed 
hyper-fine phenomena which were not well understood. Since then Oe has become possi- 
bly the most accurately determined quantity in science, being known to a precision better 
than one part per billion [2] . The corresponding theoretical calculation has achieved sim- 
ilar accuracy [3] . Because of the relatively light mass of the electron, the calculation is 
strongly dominated by QED contributions with virtual electrons, which are known to a 
good accuracy to four-loops. Using an independent determination of the fine-structure 
constant a from atomic interferometry results in a value of Ce which agrees with the 
experimental result, with an uncertainty over 30 times greater. Combining the experi- 
mental and theoretical results for Oe in terms of the fine structure constant a provides 
the most accurate available determination of a [2] . 

2 

Because of its heavier mass, —5^ — 40000, the muon anomalous magnetic moment a,, 
is far more sensitive to contributions from other sectors of the standard model, as well 
as to any potential new-physics contributions. This makes it a far more robust test of 
the standard model, and a much more interesting searching-ground for signals of new 
physics. The current experimental result, while not nearly as accurate as that for a^ is 
still remarkably precise jl]: 

a^ = 11659208.0(6.3) x 10"^°, (1.2) 

which remains a precision of better than one part per million. 

Obtaining a theoretical result for a^ of comparable precision has proved a more difficult 
task than in the case of Oe [5]. This is because, as stated above, the contributions from 
other sectors of the standard model are more significant. However the calculation has 
been brought to a point where the uncertainty is of the same order as the experimental 
uncertainty. Interestingly however, there is a discrepancy between the two values which 
exceeds the current uncertainty. This has attracted a huge amount of interest to a^ and 
lead to significant efforts to calculate contributions from potential new-physics sectors. 

The current uncertainty in a^ is strongly dominated by hadronic contributions, specif- 
ically the leading order hadronic, and hadronic light-by-light contributions. The light-by- 
light contribution has attracted significant theoretical interest, and has recently become 
the focus of considerable work using lattice simulations [HI E] . 

This work involves the leading order hadronic contribution, which we denote as alt ° , 
the best estimate of which is currently obtained by relating the hadronic vacuum polarisa- 
tion of the photon to the cross section for e~^e~ decay into hadrons, allowing a dispersive 
integral over experimental data for the cross section |8j. 



Despite the apparent accuracy of the results obtained from this procedure, there 
remain discrepancies between results from different data sets, and these discrepancies 
carry over to the total result for ajj, , and depending on the choice of method can 
result in a ~ 3cr or ~ la discrepancy between the theoretical and experimental result 
for a^. As a result, it is not clear if this method of obtaining the vacuum polarisation 
is under good control [3, [S]. In addition, the discrepancy of a^ is an important input 
for new-physics model builders when constraining their models. As such an independent 
first-principles calculation of a)j " could have a large impact on efforts to discern the 
nature of physics at and above the electro-weak scale. 

Attempts have been made to estimate this quantity using models of low energy QCD 
[H] , however for a truly robust first-principles evaluation of this quantity we must turn to 
lattice field theory techniques. 

This quantity was first tackled through lattice computation in quenched simulations 
first with domain wall fermions [10] , followed by a calculation with improved Wilson 
fermions [11]. The first dynamical simulation followed [T2| [T3] using 2-1-1 flavour staggered 
quarks, and several studies of this quantity are ongoing, using 2 flavours of improved 
Wilson fermions [Hj and twisted mass fermions [TJ]. We present a calculation oi ai 



I'fj, 
from a dynamical simulation of 2-1-1 flavour QCD with domain wall fermions. 

2 Background 

The Lande f^-factor of a fermion can be expressed in terms of the electromagnetic form 
factors Fi and F2 as 

(7 = 2[Fi(0) + F2(0)]. (2.1) 

These form factors are defined in the effective electromagnetic scattering vertex whereby 
the expression for the tree-level graph 



-^67/. (2.2) 
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p 



is replaced by its equivalent including all quantum corrections 




-ieTf,{p',p) = -ie 



l.FM') + '^m') 



(2.3) 



V p' 

From the Born approximation it can be seen that -Fi(O) = 1 to all orders, and so 

a = ^ = FM- (2.4) 



We seek to compute the effect of hadronic vacuum polarisation contributions to a^ which 
are obtained by calculating contributions to the graph in (2.3) of the form 



g, /i 




(2.5) 



{hai 



As described in [TU] the contribution to a^ from the one-loop diagram equivalent to the 



graph (2.5) with the hadronic blob removed can be expressed as 




,(1) 



a 



TC 



dQ'fiQ' 



(2.6) 



where the kernel function f{Q'^) is divergent as Q^ — ;■ and can be expressed 



f{Q')= '\ .A.v/Xa. Z{Q^) = -- V,^, ^^ . (2.7) 



l+m2Q2^(g2)2 



2mlQ^ 



From this, the expression for the hadronic vacuum polarisation contribution can be 
obtained with the insertions: 




a 



{2)had 



ay 
7tJ 



dQ'f{Q')xU{Q' 



(2.8) 



where Il{Q^) is the infra-red subtracted transverse part of the hadronic vacuum polari- 
sation 

U{Q') = n{Q') - n(0) U,,{q) = (q'g,, - q.q^Mq') (2.9) 



q, fi -^^acj-^^ (1,^ = i'n^,uiq) 



(2.10) 



at Euclidean momentum Q^ = —q^. The hadronic vacuum polarisation function n^j,(g) 
can be computed as the Fourier-transformed two-point correlator 

(2.11) 
(2.12) 



n;..(g) = j d^xe'^<--y\j,{x)Uy)) 
involving the electromagnetic current 



Jt^i^) = ^Qi'4''ii^'4'' 



where if)^ is the quark field of flavour i and Q* is its charge. The path-integral used in the 



expectation value in (2.11) will involve only hadronic fields, i.e. quarks and gluons. 



2.1 Simulation 

Our computation is performed using configurations generated by the RBC & UKQCD 
collaborations as part of their program of investigation using 2+1 flavours of domain-wall 
fermions. We investigate three lattice volumes, each with several ensembles at different 
values of the light quark mass rriu- The parameters of these ensembles are given in 
Table [TJ The ensembles at /3 = 1.75 have been generated using a dislocation suppressing 
determinant ratio (DSDR) in conjunction with the Iwasaki gauge action, with a fifth 
dimension whose extent is ^5=32 [161 E]- The lighter of these ensembles is very near to 
the physical point with a pion mass of m.,r — 180 MeV. The other ensembles used only 
the Iwasaki action and L5 = 16 
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/3 


a-i GeV 


e.n GeV^ 


aruh 


aiTLu 


243 X 64 


2.13 


1.73(2) 


0.028 


0.04 


0.02 


243 X 64 


2.13 


1.73(2) 


0.028 


0.04 


0.01 


243 X 64 


2.13 


1.73(2) 


0.028 


0.04 


0.005 


32^ X 64 


2.25 


2.28(3) 


0.05 


0.03 


0.008 


32^ X 64 


2.25 


2.28(3) 


0.05 


0.03 


0.006 


32^ X 64 


2.25 


2.28(3) 


0.05 


0.03 


0.004 


32^ X 64 


1.75 


1.375(9) 


0.018 


0.045 


0.0042 


32^ X 64 


1.75 


1.375(9) 


0.018 


0.045 


0.001 



Table 1: Parameters of the lattice ensembles used in our study. 



/5 


arriu 


Zv 


amy 


amps 


a/v 


2.13 


0.02 


0.696(2) 


0.579(6) 


0.3227(7) 




2.13 


0.01 


0.700(2) 


0.529(5) 


0.2422(5) 




2.13 


0.005 


0.699(2) 


0.505(6) 


0.1904(6) 




2.25 


0.008 


0.7380(5) 


0.388(6) 


0.1727(4) 


0.078(6) 


2.25 


0.006 


0.7385(6) 


0.366(5) 


0.1512(3) 


0.076(5) 


2.25 


0.004 


0.7387(7) 


0.356(6) 


0.1269(4) 


0.070(11) 


1.75 


0.0042 


0.664(5) 


0.570(25) 


0.1809(3) 


0.102(6) 


1.75 


0.001 


0.669(8) 


0.558(44) 


0.1249(3) 


0.105(15) 



Table 2: Relevant observables measured on our lattices. Results on the f3 = 1.75 lattices are 
preliminary and will be outlined in a forthcoming publication |17| . results for fy on the 64 x 24'^ 
lattices are currently unavailable. 



2.2 Vacuum polarisation 

We compute the lattice vacuum polarisation as 

u,Ax) = zvj2Qy{v,ix)vim, 



(2.13) 



where we have omitted the fiavour-nondiagonal terms as they contain only "disconnected" 
contributions which are expected to be sub-dominant, as will be discussed further below. 



At the sink we use the DWF conserved vector current 

^;(^) = E 2 ["^'^"^ + '^' ')(^ + i,)uixx)n^, s) - n^, s){i - ^,)u,{xmx + a, s)] 

s=l 

.(2-14) 
while at the source we have the local vector current VJ(x) = (f (x)7i,g*(x) where q^{x) = 
P^ip\x, L5 — 1) + P_ip\x, 0), and P± = |(1 ± 75). Because of the use of the local vector 
current, a factor of the vector current renormalisation constant, Zy, is included in our 
definition of the vacuum polarisation. The values of Zy used on each ensemble are given 
in Table [21 as measured in [19] . 

These correlators were generated for, and used in, the measurement of the QCD 
contribution to the electro- weak S-parameter [21]. However, they will prove perfectly 
sufficient for our purposes, as long as we are mindful of Ward Identity violations, which 
will be discussed in Sec. 12.31 

Of the two Wick-contractions arising from this correlator, we compute only the con- 
nected one. We leave the evaluation of the disconnected contribution for future work, 
but note that it is expected to be suppressed relative to the connected contribution [22] . 
This argument is also the motivation for neglecting the flavour-nondiagonal terms, and 
we will make an estimate of the systematic uncertainty that results in our conclusions. 

We Fourier transform into momentum space: 

u,M ^ZyY,Q-Y. e^''^«'(v;(x)K:(o)) (2.15) 



X 



2Txn 



using the discrete momenta g^ = ^-^ where n^ is a 4-tuple of integers, and L^ is the 
length of the lattice in the /x direction. From here, we will use the lattice momentum 

g, = ^sinf^V (2.16) 



a \ Lf, 



We associate the quantity g^ = ^„g^ with the continuum momentum Q^. 
2.3 Ward identities 



In order to ensure that this reproduces a vacuum polarisation of the form (2.9) we must 
verify that this lattice correlator satisfies the Ward identity q^H^u = which in general 
is not the case, as although both operators V and V^ have the correct continuum limit 

v;, y; ^ J^ = V^^7^# (2.17) 

the additional irrelevant operators introduced into the lattice action modify the Ward 
identity for 11^;^. In coordinate space, the Schwinger Dyson equation for II^j, reads 

where A„ is the backward lattice derivative. Because the local current used is not point- 



split, the second term in ( 2.18^ 
This is illustrated in Fig. [l 



agfi 



]) vanishes and we have as a result that e 2 g^II^j, = 0. 
where we see that it is necessary to include the factor e*^" 
in the Ward identity for the first index of II^j,, while there is no fulfilled Ward identity 
for the second index. 
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Figure 1: Illustration of Ward identity violations in II^i, on 32^ x 64 lattice at /3 = 2.25 and 
aruu = 0.004. 



2.4 Decomposing the vacuum polarisation 

We must extract from n^,y(g) the scalar vacuum polarisation n(g^) which, corresponding 



to the continuum (2.9), are related by 



^,u{q) = {q%. - q,q.Me) (2-19) 

In practice, in order to avoid any longitudinal contribution which might arise due 
to the non-conservation of Ward identities, for each momentum orientation we choose 
directions /i such that (jfj, = and compute 



U{f 



n^^(g) 



(2.20) 



where in the above there is no sum over /i. 

In Fig. |2] we show an example of the resulting vacuum polarisation function, and 
compare this to the large Q^ expansion of the three-loop continuum perturbation theory 
result from [23] , using two massless flavours of quarks and one massive flavour which we 
associate with the strange quark. This result is quoted in the MS scheme and as such we 
require the strange quark mass in our simulations expressed in MS. For this we use the 
non-perturbative renormalization factor Z^^ = 0.1533(6) (33) determined in [12]. The 
factor is quoted in the limit of vanishing light quark mass, but it is also illustrated that 
the mass dependence is extremely slight, and so we see this as satisfactory. 
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■ Lattice data 
— Perturbation Tiieory 
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Figure 2: Vacuum polarisation function J1{Q'^) as measured on 64 x 32^ lattice at /3 = 2.25 and 
arriu = 0.004. 



T^ 1 . (2)had 

Deducing a/i 



(2)had 



In order to infer the value of a/^ from our data we must carry out the integral (2.8) 
which we split into high and low momentum regions at some momentum cut Q^ 



a 



{2)had 



W 



Itl 



dQ'fiQ') X U{Q') + / dQ'f{Q') X n(Q 



(3.i: 



A continuous description of T\.{Q^) at low momenta is obtained by performing a fit 
to our lattice data, which allows us to perform the low Q^ integral. The value of n(0) 
from this fit combined with a high-momentum description of I\.{Q^) from perturbation 
theory allows us to perform the high momentum integral. As we shall see, the integral is 
strongly dominated by the low momentum contribution. 



3.1 Fitting the low Q^ region 

We have attempted to fit a continuous form to our lattice data for the vacuum polarisation 
using a number of different fit forms. The effect that the choice of fit function can have 
on the result for a/i " has been highlighted in previous studies [13], and this behoves us 
to ensure that the systematics with regard to this choice are under control. 
The suitability of a given fit-form should be judged on two main criteria: 



• Firstly, the chosen expression must describe the data closely, and must do so re- 
gardless of the range of data included in the fit. As such we require the reduced x^ 
of the fit to be consistently low as a function of Q^ which defines the range of data 
in the fit. 

• Secondly, in order to deduce that the fit-form results in an integral over momentum 
which is relatively stable, we desire that the result for a/^ is again relatively 
stable as a function of Qq. 

Ref. jT2| also illustrated the use of a fit form originating in the expression for the 
vacuum polarisation calculated in chiral perturbation theory. The dominant component 
of this expression is due to the vector meson contribution, which at tree-level is 



2_g 
SQ-' + ml 



^'rm = -.V^^ (3-2) 



where the vector decay constant fy is defined 

(Q|J^|V,p,e) =my/ve^(p). (3.3) 

Motivated by this expression the fit-form we use is closely related, differing only in 
the inclusion of the contribution of an additional vector resonance, 



^m = A- —-^ - — ^^. (3.4) 



Fl Fl 

Q'^ + ml Q^ + ml 

The one-loop contribution from the pseudoscalar sector, shown in [13] to have small 
momentum dependence, will not strongly affect our results and so, in our effort to make 
a continuous description of the lattice data, it will be omitted from our fit ansatz. 
We fit the lattice vacuum-polarisation data in two ways: 

• Firstly using A, Fi^2 and ?7ii^2 as free parameters. 

• Also, fixing the parameter nii to the mass of the vector meson my as measured in 
J19j . This we do by constraining mi to lie in the one-sigma band defined by the 
estimate of ttiy and its variance. This method was found to maintain the stability 
of the fit routine, while incorporating the extra information provided by my. In 
this fit y4, Fi^2 and m2 remain as true free parameters. 

The behaviour of such fits are shown in Fig. [3j Clearly such a form is a very good 
representation of the data, over practically the whole range of Q}j. In addition the 
results for a/^ "" using such fits are very stable as the fit range is varied, allowing far 
greater confidence in the reliability of the result. In particular we conclude that using 



a fit form (3.4) with the mass of the first pole fixed to the ground-state vector meson 
mass to be the optimal method of describing the lattice data for the hadronic vacuum 
polarisation. 



In Fig. |4]we see the value of the fit parameter mi from (3.4) as determined from fits to 



the lattice vacuum polarisation. The value of my obtained in [19j is shown in green, and 
this defines the band in which mi was constrained to reside in the fixed version of this fit. 
We have not attempted to model 0(4) breaking effects present in our data. Though such 
effects do appear to be present to a moderate degree on certain ensembles, they do not 
prevent the extraction of a reasonable signal from our data at this point. These effects 
could also be alleviated by the use of twisted boundary conditions [24j. 
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Figure 3: Properties of fits to the lattice vacuum polarisation using the ansatz (3.4) on the 
P = 2.25 lattice at arriu = 0.004. 



3.2 Evaluation of (3.1) 



Illustrations of the integrand can be seen in Fig. [5j Because the integrand is dominated 
by contributions in the low momentum region, we change our integration measure to 
better sample the region of interest. To do this, we make the change of variables 



t 



1 



1 + log^ 
and so the integral over the low-momentum region becomes 



(3.5) 



dQ^fiQ^) X fi(g' 



dtfiQ^)xn{Q' 



X 



t2 



(3.6) 



Overlaid on the depiction of the integrand in Fig. [5] is the appropriately subtracted 
and rescaled vacuum polarisation data. We see from this that, while a large portion of 
the constraint on the fit is consistently derived from data at higher momentum, the fit is 
always consistent with the data at low momentum, the region where the integral receives 
the dominant contribution. 

In particular in Fig. 5b we see that on the larger lattices at /3 = 1.75 using the 
Iwasaki+DSDR action, the data point at the lowest momentum sits exactly where the 
integrand reaches a maximum, and there are numerous data points in the dominant 
region, constraining the fit. Clearly using lattices of such size will help in obtaining 
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Figure 4: Value of the fit parameter ami in fits using the ansatz (3.4) on the /3 
at arriu 



2.25 lattice 



0.004. The vector mass arriY as determined on this lattice is shown in green. Note 
in the fit where mi was fixed, it was only constrained to lie within the green band. It is clear 
that for a high Q^, ttt-i will emerge at the upper limit of the band, indicating some tension 
between the fit-form and the data, but as can be seen in Fig. [3j this has very little impact on 
the goodness of the fit. 



a precise result for this quantity, and this must be combined with the use of twisted 
boundary conditions [H] in order to access data at lower values of the lattice momentum. 




1/(1+log(Q '/q')) 
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Figure 5: Examples of the integrand in the rescaled integral (3.6). 
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4 Results 



We extract our final results from the fit using (3.4) with the first mass fixed to that of 



the vector meson as measured on each ensemble. Observing the behaviour of the reduced 
X^ as the fit range is varied, we choose a suitable value for Q^ for each ensemble which 
provides the most reliable result. We attempt to choose a cut which provides a low 
reduced x^ preferably where the parameter mi agrees without tension with my. This 
produces the results shown in Table JSJ where we also quote the reduced x^ of the fit, and 
the resulting values of the remaining associated free parameters. 

These results are also shown as a function of m^ in Fig. p[ where we compare them 
to previous 2+1 flavour results from [13]. Also shown is an extrapolation to the physical 
point, using a quadratic chiral ansatz. This produces a final result for the leading order 
hadronic vacuum polarisation contribution the anomalous magnetic moment of the muon 
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Figure 6: Integrated result for a^ as a function of the pseudoscalar mass squared. 
We have also investigated the effect of modifying the kernel function in the integrand 



(2.8) in the manner outlined in [15], where in an effort to moderate the variation of the 



outcome of the integral as a function of the quark mass, the momentum argument of the 
kernel function is rescaled by a factor of the ratio of the value of a relevant observable H 
(the mass of the vector meson appears to be an optimal choice) measured at the simulated 
quark mass to its physical value. This effectively defines the calculation of a new quantity 
which approaches the desired a)j in the physical limit. We show the results of such 
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a calculation in Fig. 7a, along with an accompanying cliiral extrapolation. The chiral 
variation in this redefined quantity is such that it allows for a linear extrapolation in 
quark mass. For the lightest point in our simulation we include the unmodified result 
outlined in Table |3] since for this ensemble the measured vector mass mv is consistent 
with the physical value. This method does indeed moderate the chiral behaviour of the 
result, however it has little effect on our data at light quark masses, primarily because 
the lattice vector meson masses are very near that of the physical p meson, and, as of 
now, are not determined to any great precision on these lattices. As such this technique 
does not improve our chiral fit at this time, producing a compatible result with a similar 
uncertainty: 



a 



{2)had 



605(24) X 10 



-10 



(4.2) 



In Fig. 7b we compare both chiral extrapolations, with H = 1 denoting the standard 



method, and H = my indicating the modified prescription of [15] using the vector mass 
my. 





(a) Results using modified prescription H = my- (b) Comparison of results from both methods. 

Figure 7: Analysis of results for modified prescription using H = my- 



In Fig. [sl our result ( |4.1 ) including a 5% statistical error arising from the chiral 
extrapolation (estimated from the discrepancy between (4.1) and (4.2) ) is compared to 
recent 2+1 fiavour lattice results [13jj along with recent result arising from dispersion 
integrals over experimental data from scattering data. Our result is in rough agreement 
with other results bearing in mind that we have neglected the disconnected contributions 
to our correlators, producing a systematic deviation which is bounded to be of the order of 
10%. At this time we are not in a position to improve on previous theoretical evaluations 
oi ttfi " and so cannot comment on the scale of the discrepancy between the measured 
value of a^ and the Standard Model prediction. However it is clear that the next iteration 
of this calculation with planned improvements is very likely to be in a position to begin 
clarifying the situation concerning this discrepancy. 



In Tabled we attempt a comparison of the value of Fi (defined in (3.4) ) resulting 
from our fit, to the vector decay constant as measured on each lattice, according to the 
relation expressed in (3.2). Note, we do not have a result for fv on the 64 x 24^ lattices 



at this time, although the ratio of the vector coupling to the vector and tensor currents 
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Figure 8: In Fig. 



we compare recent results for a}j . 



Our result is compared to previous 
2+ If lattice result [13j and results from a dispersive integral over cross-section data from e~'"e~ 
data [H] and from e+e— and r data [23]. In Fig. 8b we compare the full result for a^ including 



the results corresponding to the dispersive results in Fig. [8a] and the current experimental result 
from [41. The increment between the labelled ticks on the horizontal axis is 10~^ 



was studied in 



We also make the comparison suggested by the one-loop correction 



where 



to (3.2) as computed in [13] whereby the relation F^ ~ |/y is replaced by F^ 



/. 



C^ 



6 



(4vr/.: 



m^log 



mz 



+ m^ 



log 



ni 



K 



2x^2 



(4.3) 



/^ / V /^ 

with 171^^ and uik the pion and kaon meson masses, f^^ the pion decay constant, and /i the 
chiral scale, taken as 1 GeV. We note that in our fits we have not included the one-loop 
contribution from the pseudoscalar sector, and so this comparison can only serve as a 
rough indication, and we do not necessarily expect close agreement. We observe that the 
value of Fi emerging from our fits is on the correct scale when compared to the measured 
/y, simply supporting the credibility of our fits. 



5 Conclusions 

We present a fully dynamical calculation of the leading-order hadronic vacuum polarisa- 
tion contribution to the anomalous magnetic moment of the muon, using a 2+1 flavour 
simulation lattice QCD using domain wall fermions. Although we have an expensive 
fermion discretisation, we improve the accuracy of our result by convolving an accurate 
determination of the ground-state vector meson mass with our determination of the lattice 
hadronic vacuum polarisation in order to suppress the systematic uncertainty associated 
with the choice of fit ansatz. Our chiral extrapolation involves lattices at different bare 
couplings, and thus different lattice spacings, however at this level of precision we do not 
detect any significant discretisation, or finite volume errors in our result. Our final result 
we take to be 

^(2)had ^ 641(33)(32) X 10"^° (5.1) 
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where the first error is statistical and the second is an estimate of the systematic error 
arising from the extrapolation to the chiral limit, taken as 5%, motivated by the varia- 
tion between the results (4.1) and (4.2). Our largest systematic uncertainty arises from 
the omission of the disconnected contributions and is of the order of 10% [27] . In order 
to obtain a more comprehensive and accurate result, we must include the disconnected 
contributions in our calculation. Furthermore, this being a first effort at deducing this 
quantity from our lattices, we have plans to improve it in a number of ways. In addi- 
tion to the enhancement of our statistics, we would like to obtain a higher momentum 
resolution through the use of twisted boundary conditions, and also to explore the use of 
stochastic sources to further enhance our signal. With these improvements we would hope 
to decrease the uncertainty in our result significantly and thus begin to clarify the true 
discrepancy, if any, between the actual value of a^ and its prediction from the Standard 
Model. 



/3 


arriu 


Ql GeV2 


n.d.f 


a}f^ X 10^0 


aFi 


ar7i2 


aF2 


2.13 


0.02 


4 


0.38(17) 


345(16) 


0.114(4) 


1.48(19) 


0.31(5) 


2.13 


0.01 


3.5 


0.07(6) 


430(22) 


0.110(4) 


1.50(23) 


0.32(7) 


2.13 


0.005 


3.5 


0.14(5) 


436(50) 


0.097(14) 


1.16(18) 


0.24(3) 


2.25 


0.008 


6 


0.18(11) 


452(23) 


0.079(2) 


1.14(4) 


0.26(1) 


2.25 


0.006 


6 


0.10(6) 


484(33) 


0.075(3) 


1.07(7) 


0.24(2) 


2.25 


0.004 


9 


0.06(3) 


568(29) 


0.079(2) 


1.23(3) 


0.28(6) 


1.75 


0.0042 


2.5 


0.16(9) 


536(36) 


0.108(20) 


1.27(20) 


0.26(3) 


1.75 


0.001 


2.5 


0.27(13) 


646(55) 


1.06(11) 


1.58(61) 


0.37(27) 



Table 3: Results for the hadronic contribution to the muon anomalous magnetic moment. 



p 


arUu 


fv MeV 


^Fi MeV 


yifMeV 


2.13 


0.02 




242(10) 


179(7) 


2.13 


0.01 




234(8) 


166(6) 


2.13 


0.005 




205(30) 


144(20) 


2.25 


0.008 


178(13) 


221(6) 


155(5) 


2.25 


0.006 


174(11) 


211(10) 


147(7) 


2.25 


0.004 


160(26) 


222(5) 


155(4) 


1.75 


0.0042 


140(9) 


192(27) 


129(19) 


1.75 


0.001 


144(20) 


179(18) 


127(12) 



Table 4: Comparison of the vector decay constant as measured on our lattices, to the amplitude 
of the lowest resonance contribution emerging from our fit to the lattice vacuum polarisation. 
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